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ABSTRACT 

We perform a simple test showing that accreting halos grow inside-out and have density pro- 
files that are statistically indistinguishable from those of halos having undergone major mergers. 
The test also reveals that the ability for the NFW analytical expression to fit the density profile 
of halos spoils towards large masses at any redshift. On view of these results, we analyze the 
M-c relations on extreme mass and redshift domains predicted by usual toy models based on 
numerical data and one physical model based on the two assumptions above confirmed by the 
test. 

Subject headings: cosmology: theory — dark matter — galaxies: halos 



1. INTRODUCTION 

A^-body simulations show that dark matter ha- 
los have spherically averaged density profiles that 
can be fitted by the NFW law (Navarro et al. f 997) 
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where Tg and ps are the halo scale radius and char- 
acteristic density, respectively. The Einasto pro- 
file (Navarro et al. 2004; Merritt ct al. 2005, 2006; 
Navarro et al. 2010) has also been shown to give 
similar or even better fits to halo density profiles. 
But the correlations between the NFW shape pa- 
rameters, namely between the halo mass M and 
the concentration c, defined as the halo radius R 
over Tg, or between rg and the mass Mg = M{rs) 
encompassed by it, at any given redshift (z) still 
play a crucial role in the modeling of structure 
formation. 

Much effort has been done in the last decade 
in accurately determining these relations from TV- 
body simulations (e.g. Zhao et al. 2009; Klypin 
et al. 2010; Muiioz-Cuartas et al. 2010; Prada et 
al. 2011 and references therein). However, such 
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numerical relations are not very practical because 
they refer to only a few redshifts and cosmolo- 
gies. Moreover, they face the fundamental prob- 
lem that numerical simulations have a limited dy- 
namic range, which severely restricts the mass 
and redshift domains where they can be inferred. 
To circumvent this problem, toy models are often 
built that fit the numerical data and are supposed 
to give acceptable predictions beyond the domain 
covered by simulations. 

On the other hand, the NFW law does not ex- 
actly fit the halo density profile, which may in- 
troduce poorly known biases in the best-fitting 
Tg and Ps values depending on the specific halo 
sample definition and fitting procedure (including 
the radial extent) used. Tasitsiomi et al. (2004) 
found indeed that very massive halos at any given 
z have density profiles that become closer to a pure 
power-law with minus index equal to ^ 1.5-2.0. 
According to these authors, this would be the con- 
sequence that massive halos aggregate mass at a 
lower rate, through minor mergers contributing to 
smooth accretion, than less massive halos, which 
suffer major mergers more frequently. But this in- 
terpretation has not been proved. The idea that 
the density profile of halos depends on the rate 
at which they aggregate mass is often invoked, in- 
deed, to explain some of the observed trends of 
halo density profiles (e.g. Wechsler et al. 2002; 
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Zhao et al. 2003; Munoz-Cuartas et al. 2010). 
However, other results of numerical simulations 
show the opposite result that there is no particu- 
lar signature of major mergers in the halo density 
profiles (e.g. Huss, Jain & Steinmetz 1999; Moore 
et al. 1999; Wechsler et al. 2002; Nusser & Sheth 
1999; Romano-Diaz et al. 2006). 

Clearly, the existence of a reliable physical 
model of halo structure could shed light on this 
apparent conflict and might be used to predict ac- 
curate M-c or Tg-Ms relations at any M and z 
domain as well. Salvador-Sole et al. (2007; here- 
after SMGH) built one simple model that relies on 
the assumption that halos develop from the inside 
out during accretion phases. The typical density 
profile for purely accreting halos is therefore be set 
by the typical accretion rate provided by the ex- 
cursion set formalism, which leads to density pro- 
files a la NFW except at very large masses where 
they tend to deviate from the NFW profile in the 
direction found by Tasitsiomi et al. (2004). 

The inside-out growth of accreting halos is 
supported by the results of iV-body simulations, 
which show that halos evolve by conserving the 
radial mapping of about 80 % of their particles 
(Wang et al. 2010 and references therein). The 
fact that 20 % of them do not preserve their ra- 
dial mapping might simply be due to the effects of 
major mergers yielding the full rearrangement of 
halos. But this is hard to prove because of the dif- 
ficulty to restrict the analysis to purely accreting 
haloqj. But even if this were the case, it would not 
be clear why the SMGH model is able to predict 
the right NFW profile of halos without including 
major mergers. 

Recently, Salvador-Sole et al. (2011a) have 
shown from theoretical arguments that accret- 
ing halos must grow, indeed, from the inside out 
and that the density profile of purely accreting 
halos must be indistinguishable from that of halos 
having suffered major mergers. Hence, the typical 
halo density profile predicted by the SMGH model 
is fully justified despite the real growth of these 
objects both through accretion and major merg- 
ers. Unfortunately, those theoretical arguments 
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are not straightforward, so the checking against 
A'^-body simulations of these important results is 
necessary. 

In the present Letter, we perform one simple 
test using numerical data which confirms that: i) 
halos grow inside-out during accretion phases, ii) 
the density profile of purely accreting halos is in- 
distinguishable from that of halos having suffered 
major mergers and iii) the fit by the NFW law of 
halo density profiles depends on Ad (at any given 
z), which causes some trends often misinterpreted 
as proving that the profiles of halos depend on 
whether they suffer major mergers or nolH.The 
implications of these results on the M-c relations 
predicted far away from the domain covered by 
numerical data is then analyzed. 

2. THE TEST 

If halos had density profiles strictly of the NFW 
form and grew inside-out by pure accretion, the 
values of the rg and Mg adjusting their density 
profiles at different epochs would clearly not vary 
and the rs-Ms relation itself would remain un- 
changed. Actually, halos rearrange in major merg- 
ers, so the and Ms values suddenly change in 
those events. But, if the density profiles emerging 
from major mergers were indistinguishable from 
those of purely accreting halos, the whole r^-Ms 
relation should be kept unaltered regardless of the 
aggregation history of halos. That is, it should 
be time-invariant. Moreover, as the points of the 
r^-Ms relation would not depend on halo mass, 
the Tg-Mg curve should not privilege any partic- 
ular mass. Hence, the rg-Mg relation should also 
be scale-free. 

Conversely, if accreting halos did not have den- 
sity profiles of the NFW form or they did not grow 
inside-out or still the profile for purely accreting 
halos were distinct from that of halos suffering ma- 
jor mergers, the Tg and Mg values would change 
as they evolve and the r^-Mg relation would vary 
with z. Moreover, as major mergers dominate 
the growth of low-mass halos and smooth accre- 
tion that of very massive ones, the rg-Mg relation 
would not be scale-free. 



•^This is not to be confused with the fact that the density pro- 
file of individual purely accreting halos depends, of course, 
on their past accretion rate, used in some halo age esti- 
mates. 
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Fig. 1. — Typical rg-Mg relations at z = 0, 
1, 2 and 3 according to Zhao et al. (2009) (red 
lines), Klypin et al. (2010) (blue lines) and Munoz- 
Cuadras et al. (2010) (green lines) and predicted 
by the SMGH model (black hnes). The straight 
dashed black lines trace the ideal rg-Mg relation 
defined in section |3l 

Thus, the time-invariant, scale-free behavior of 
the Tg-Ms relation is an unambiguous proof of the 
inside-out growth of accreting halos and the sim- 
ilarity of density profiles for purely accreting and 
merging halos. Hence, it is an unambiguous proof 
of the validity of the SMGH model (as far as the 
halo mass growth is reasonably well-described by 
the excursion set formalism). 

In Figure [1] we show the rg-Mg relations 
drawn from the toy models of M-c relations (see 
Sec. [21) proposed by Zhao et al. (2009; top-right 
panel), Klypin et al. (2010; bottom-left panel) 
and Muhoz-Cuartas et al. (2010; bottom right- 
panel) fitting numerical data in their simula- 
tions of the same concordance cosmology with 
inm,nA,h,as) = (0.27,0.73,0.70,0.82). (Muhoz- 
Cuartas et al. actually use a slightly different 
concordance model, but this makes no difference 
at the resolution of the Figure.) The curves at 
z = are very close to straight lines in log-log as 
expected, but, for increasing z, they progressively 
bend upwards, so they do not overlap. 



Does this mean that the SMGH model is 
wrong? Not really. In the same figure (top-left 
panel), we plot the Vg-Mg curves predicted from 
the SMGH model. Again, the log-log relation at 
z = is very approximately a straight line, but, 
as z increases, the relations progressively bend 
upwards and the curves do not overlap. As in the 
SMGH model all halos accrete (there is no major 
merger) and develop, by construction, from the 
inside out, this can only be due to the fact that 
they do not strictly have density profiles of the 
NFW form. Indeed, even if the density profile of a 
halo does not vary during accretion, its fit by the 
NFW law will then result in best values of Vg and 
Ms that slightly shift as the radial extent of the 
halo (setting the fitting radial range) increases. 

This interpretation of the actual behavior of the 
rg-Mg relation is confirmed by the results plotted 
in Figure [2] In panel (a) , we show again the rg- 
Mg relations predicted at several z's by the SMGH 
model but in a somewhat smaller rg range and in 
panel b) we show the tracks followed for decreas- 
ing z by a few individual halos tracing such Vg- 
Mg relations. In the SMGH model, all halos de- 
velop inside-out. Yet, the halos still move in the 
Tg-Mg plane along one universal direction due, as 
checked, to the effect just mentioned. Remarkably, 
the slope of that universal shift is identical to that 
{y = a;3xO-55') Qf ^]^g shift also undergone by halos 
in the simulations by Zhao et al. (2009) (see their 
Figure 22). Moreover, the shifts found by Zhao et 
al. end up at z = along the same straight line 
tracing the SMGH rg-Mg relation at the same z. 
This unambiguously proves that such shifts have 
the same origin in the SMGH model as in Zhao et 
al. simulations. Zhao et al. interpreted them as 
due to a progressive change in the inner halo den- 
sity profile during slow accretion phases (see also 
Zhao et al. 2003). However, as revealed by the 
SMGH model, the inner density profile for these 
halos does not actually change; the small shift in 
the rg-Mg plane is due to the varying best-fitting 
rg and Mg values obtained as the halo radial ex- 
tent increases due to the fact that at massive halos 
are not well-fitted by the NFW profile. 

When this effect is corrected for, that is, when 
the values of the shape parameters adjusting the 
density profile at any z are replaced by those ad- 
justing the same profile extended so to encompass 
the whole radial range at z = 0, then the log-log 
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Fig. 2. — Panel (a): Vg-Ms relation (full lines) 
predicted by the SMGH model at z = 0, 2, 4, 
6 and 8 (respectively in black, green, red, blue 
and magenta). Panel (b): Tracks followed by typi- 
cal accreting hales (dashed red lines) with current 
masses equal to 10", lO^^, IQi^, lO^^ and 10^^ Mq 
(from bottom to top) and the same redshifts (cir- 
cles) as in the left panel. In red dashed straight 
line, the universal direction followed by slowly ac- 
creting halos found by Zhao et al. (2009). In full 
black line, the a straight line with the same slope 
as the Tg-Mg relation predicted by the SMGH 
model at z = 0. 



rg-Mg relations at any z become straight lines that 
overlap. In other words, the test confirms that ac- 
creting halos evolve inside-out and that the den- 
sity profiles of purely accreting halos are statisti- 
cally indistinguishable from those of halos having 
suffered major mergers. 

3. IMPLICATIONS FOR THE M-c RE- 
LATION 

Integrating 47rr^ times the NFW profile (eq. [1]) 
written in terms of Ms instead of ps and taking 
into account the relation Ms(rs) = Ar^ with A = 
9.99 X lO^"* M0 Mpc""" and v = 2.481 fitting the 
numerical r^-A'U relations at z = 0, one is led to 
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where Avir(z) and p{z) are respectively the halo 
overdensity and mean cosmic density at z. Equa- 



tion ([2]) therefore defines the M-c relation for ha- 
los strictly endowed with the NFW profile (in par- 
ticular at z = 0), growing inside-out during accre- 
tion periods and with no signature of major merg- 
ers. Such an M-c relation is from now on called 
"the ideal M-c relation" . 

In Figure [U we compare such an ideal M- 
c relation with the M-c relations predicted by 
the SMGH model and the toy models by Zhao 
et al. (2009), Klypin et al. (2010) and Mufioz- 
Cuartas et al. (2010) derived from numerical data, 
respectively given by 
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with t{z) equal to the cosmic time corresponding 
to z and to.oiiM) the time when the main halo 
progenitor gained 4% of its mass M at z. 
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where 5{z) is the critical overdensity for collapse at 
z, a(z) = 0.029Z-0.097 and b{z) = -110.001/(z-H 
16.885) + 2469.72/(z + 16.885)^. 

All the curves drawn from numerical data be- 
have in a similar way. Contrarily to the rg-Mg 
relation, they show a strong dependence on z. At 
z = they closely follow the ideal M-c relation 
but, at larger z, they progressively deviate from 
it at large masses likely due to the effect men- 
tioned in section [2] In fact, for large enough M, 
the ideal decreasing trend brakes and the empiri- 
cal curves level off or even begin to increase again 
(in the case of Klypin et al. 2010; see also Prada 
et al. 2011). The M-c relation predicted by the 
SMGH model shows similar trends although it re- 
mains closer to the ideal M-c relation until slightly 
larger M at each z. This might refiect a slight 
flaw of the SMGH model at very large M. But 
the possibility cannot be ruled out that the em- 
pirical M-c relations are somewhat biased owing 
e.g. to the impact of using non-fully relaxed halos 
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Fig. 3. — Comparison between the M-c rela- 
tions predicted by the SMGH model (full black 
lines) and proposed by Zhao et al. (2009) (red 
lines), Klypin et al. (2010) (blue lines) and Muiioz- 
Cuadras et al. (2010) (green lines) dX z — Q, 1, 
2 and 3 in the same cosmology as in Figure [T] 
We distinguish that part of the curves covered by 
simulations (thick colored lines) and their extrap- 
olations outside it (thin colored lines). The ideal 
NFW M-c relation (see text) is in dashed black 
lines. 

(more abundant at large masses) or of estimat- 
ing the typical halo concentration in mass bins by 
means of the median value. 

In Figure SI we plot the same M-c relations 
extended over a much wider domain (M down 
to 1 Mq and z up to 15). We exclude the M- 
c relations by Muhoz-Cuartas et al. (2010) be- 
cause these authors state explicitly that their toy 
model is a fitting expression that cannot be ex- 
trapolated. As can be seen, the two sets of ex- 
trapolated empirical M-c relations have a lower 
bound at about log{c) ^ 0.6 contrarily to the ideal 
M-c relations. However, in Klypin et al. (2010), 
the curves increase again after reaching that min- 
imum, whereas, in Zhao et al. (2009), they level 
off. Notice that this discrepancy affects the part 
of the curves that adjusts numerical data, mean- 
ing that the disagreement between both groups is 
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Fig. 4. — Same as Figure [3] (and identical symbols) 
but for a much wider mass and redshift range. The 
different curves in each panel are for z equal to 0, 
1, 3, 5, 9, 15 (from top to bottom on the left of 
each panel) . 

very relevant. At small Af , the two sets of curves 
diverge from each other even more markedly: in 
Zhao et al. they show flat asymptotes, while in 
Klypin et al. they have increasing asymptotes sim- 
ilar to the ideal M-c relations, although their de- 
pendence with z is substantially less marked. 

The z-dependence of the SMGH curves is in- 
stead quite similar (for intermediate M at each z) 
to that of the ideal M-c curves. Moreover, the M- 
dependence at any z is also much closer to that of 
the ideal M-c relations than for the extrapolated 
empirical relations. They only substantially de- 
viate from each other at very small M . In brief, 
the SMGH M-c relations are globally in better 
agreement with the ideal M-c relations than the 
pseudo-empirical M-c relations. 

Certainly, the ideal relations must not necessar- 
ily trace the M-c relations of real halos because, as 
shown, their density profiles are not exactly fitted 
by the NFW law. This translates into an appar- 
ent variation of rg and c with increasing z, which 
is not present in the ideal relations. Thus, on the 
sole basis of the previous comparison, it is hard to 
conclude which M-c relation is the most reliable 
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for real simulated halos in such extreme M and z 
domains. As mentioned, the SMGH model might 
be affected by some slight flaw at the large mass 
end (and perhaps at the small mass end as well) 
at any given z. But the empirical curves might be 
biased too, as suggested by their distinct behav- 
ior at very large masses where numerical data are 
available. More importantly, their extrapolation 
at very small M show so different behaviors at any 
z that they cannot be trusted. Last but not least, 
Salvador-Sole et al. (2011b) have recently shown 
that the properties of halo substructure found in 
high-dynamic range TV-body simulations (Springel 
et al. 2008), very sensitive to the behavior of the 
M-c relation in an M and z domain as wide as 
the one considered here, clearly favor the SMGH 
model in front of the toy models by Klypin at al. 
and Zhao et al. 

4. SUMMARY 

We have performed a simple test showing un- 
ambiguously that accreting halos grow inside-out 
and that the density profiles of halos grown by 
pure accretion are indistinguishable from those of 
halos having suffered major mergers. The varia- 
tion in the NFW shape parameters as halos ac- 
crete is due to the fact that the density profile for 
very massive halos is not well-fitted by the NFW 
law, which causes the best-fitting values of rg and 
Ps (or Ms) to slightly vary as the radial extend of 
halos increases. After correction of this effect, the 
Ts-Ms relation is linear in log-log and overlaps for 
all z as expected. 

Over very wide M and z domains, the ideal 
M-c relations for halos with exact NFW profiles 
and evolving inside-out by accretion deviate no- 
tably both from the extrapolated M-c relations 
inferred from numerical data as well as from the 
M-c relation predicted by the SMGH model. This 
difference is once again partly due to the fact that 
the density profile for simulated and modeled halos 
is strictly not of the NFW form. However, the M- 
c curves predicted by the SMGH model stay closer 
to such ideal M-c curves over a wider M and z do- 
main than the extrapolated M-c relations drawn 
from toy models, which also notably deviate from 
each other. This suggests the better behavior of 
the SMGH M-c relations at such extreme M and 
z domains, a conclusion that meets that found on 



the basis of completely independent arguments by 
Salvador-Sole et al. (2011b). 
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